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Abstract. We prove a regularity result for the Poisson problem -An = /, 
"lap = g on a polyhedral domain P C M 3 using the Babuska-Kondratiev spaces 
/C™(P). These are weighted Sobolev spaces in which the weight is given by 
the distance to the set of edges El 1331 . In particular, we show that there is 
no loss of /C^-regularity for solutions of strongly elliptic systems with smooth 
coefficients. We also establish a "trace theorem" for the restriction to the 
boundary of the functions in /C™(P). 
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1. Introduction 

Let tt C R n be a smooth, bounded domain. Then it is well known P ITfil YM 1471 
153) that the equation 

(1) Au = f € H™- 1 ^), ii=0 on 9fl, 

has a unique solution u £ H m+1 (£l). In particular, u will be smooth on f2 if 
/ is smooth on fL This well-posedness result is especially useful in practice for 
the numerical approximation of the solution u of Equation l(T|. see for example 
P 1121 ITfi| among many possible references. 

In practice, however, it is rarely the case that Q, is smooth. In fact, if dft is 
not smooth, then the smoothness of / on f2 does not imply that the solution u 
of Equation Q is also smooth on fl. Therefore there is a loss of regularity for 
elliptic problems on non-smooth domains. Wahlbin (see also 0EH1IEE!) has 
shown that this leads to some inconvenience in numerical applications, namely that 
a quasi-uniform sequence of triangulations on f2 will not lead to optimal rates of 
convergence for the Galerkin approximations Uh of the solution of Q . 
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The loss of regularity can be avoided, however, if one removes the singular points 
by "sending them to infinity" by suitably changing the metric with a conformal 
factor. It can be proved then that the resulting Sobolev spaces are the "Sobolev 
spaces with weights" considered for instance in 111! 1121 l3~3"| and in several other 
papers. A related construction, leading however to countably normed spaces, was 
considered in |29| . Let / > be a smooth function on a domain f2, then the mth 
Sobolev space with weight f is defined by 

(2) *™(fi; /) := {u, f^~ a d a u G L 2 (n), \a\ < m}, m G Z+, n£l 

The regularity result for Equation 0} extends to polyhedral domains P in three di- 
mensions with the usual Sobolev spaces replaced by the spaces /C" l (P) := /C™(P; <&), 
i? being the distance to the edges. The spaces /C™(<9P;$) on the boundary are 
defined similarly for m G Z + := {0,1,...}; for m G M+ they are defined using 
interpolation. 

Theorem 1.1. Let P C M be a polyhedral domain. Let m G Z + and a G R. 
Assume that u G /C^ +1 (P), Au G ^"^(P), and u| SP G JC"'^ {dP , d) , then u G 
^a+i (^) a7l( ^ i/iere exists C > independent of u such that 

IMI^YYtr) ^ c '(II Au IIk™_- 1 (p) + IMIk° +1 (p) + IKIspll^+v^ar,*))- 
TTie same result holds if we replace A wii/i a strongly elliptic operator or system. 

Theorem 1 1.1 1 is well known in two dimensions, i. e., for polygonal domains, and 
for domains with conical points |15l l3"3") . See also |231 1341 13"5] where similar results 
were proved using a dyadic partition of unity technique. For the result in two 
dimensions, d is the distance to the vertices of the polygonal domain considered or 
to the conical points. In general, in d dimensions, one takes t?(x) to be the distance 
to the set of non-smooth boundary points of P. Significantly less papers have dealt 
with the case of three dimensions. Nevertheless, let us mention the following. A 
general and far reaching theory (valid also in higher dimensions) was developed by 
Dauge in |24] . Regularity estimates based on singular function expansions were 
proved by Apel and Nicaise and Lubuma and Nicaise [HSj. These results were 
then applied in these papers in order to obtain optimal rates of convergence in the 
Finite Element Method. In 39 , Mazya and Rossmann have obtained similar results 
using estimates on Green functions. Buffa, Costabel, and Dauge ^H] have proved 
or stated similar regularity and well posedness results for polyhedral dimensions in 
three dimensions. Our modified weight rn was introduced in |22| . where the above 
regularity theorem was proved for m = 1. In |32| . Kellogg and Osborn have obtained 
regularity results of a similar kind for the Stokes operator. Borsuk and Kondratiev 
established many regularity results for Dini-Liapunov regions in K", n > 3, in 
their recent monograph 37 . Note that the notion of a Dini-Liaponov region is a 
generalisation of a domain with C 1 ' Q -bondary. See also pfl I2T1 1231 ETT1 ETT11 \M US] . 
to mention just a few other papers. A regularity result valid in all dimensions was 
obtained in P using "Lie manifolds." 

We are grateful to one of the referees, who pointed out to us that Theorem ll.il 
can also be obtained from the results of the monograph [46j. In this paper, we 
follow pQ, but we use more elementary methods that lead to a short proof. We also 
introduce some ideas that are specific to polyhedral domains in three dimensions 
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and may be useful in applications to Numerical Analysis. Moreover, our paper is 
self-contained and the references to are only for comparison. 

We would like to stress that Theorem 11.11 does not constitute a Fredholm (or 
"normal solvability") result, because the inclusion JCJJS (P) — > ^a+i(P) ^ s n °t 
compact for all m and a pQ. By contrast, if P is a polygon, then P = —A with 
Dirichlet boundary conditions is a Fredholm operator from IC^X (P) to /C^I^ 1 (P) 
precisely when a is different from kn/a, where k £ 7*. k =/= 0, and a ranges through 
the angles of the polygon [331 134) . 

The Poincare inequality ||u||x;i(p) < C f ||V , u|| i 2(Q) proved in [T3|, gives that A is 
coercive on the space JC\(V) and hence the map A : /C*(P) tl {u = on 3P} — > 
/Cl^(P) is a continuous bijection. By combining this with Theorem ll.il we obtain 
that 

(3) A : /C™Y (P) n{^Oon 5P} -> /C^ 1 (P) 

is a continuous bijection, for any m G Z + and |a| < 77, with 77 depending only on P. 
The same result holds if A is replaced with P + cp, where P is a strongly elliptic 
system with smooth coefficients and cp > and 77 > are constants depending 
only on P 14 . 

To prove Theorem ll.il we first introduce the weighted Sobolev spaces /C™ (dF, $) 
on the boundary of P. For m Z + , these spaces are defined by duality and interpo- 
lation. Then we provide an alternative definition of the spaces /C™(P) := /C™(P, 1?) 
and IC™(dP,'d) using partitions of unity. This allows us to define a trace map 
/C™(P) -> /C^_" 1 1 / 2 2 (aP,'i?), which extends the restriction map and is a continuous 
surjection, as in the case of a smooth domain. We also show that any differ- 
ential operator P of order m with smooth coefficients induces a continuous map 
P : /C*(P) - /C'I£(P). 

We need to introduce an enhanced space of smooth, bounded functions C°°(EP), 
which contains the cylindrical and spherical coordinates functions and is minimal 
with this property. In particular, C°°(n) C C°°(SP) C C°°(fi). Let p P (p) be the 
distance from p to the vertex P of P and r e (p) be the distance from p to the line 
determined by the edge e of P (for P non-convex we need to slightly change the 
definition of r e ). Then pp,p e £ C°°(SP), although they are not smooth functions 
on fl in the usual sense. Let A and B be the end vertices of the edge e (i. c, 
e = [AB]). We further define r e :— p~j^p~^r e and rp = Y\ e r e x Jlp/ ^- Then 
f e ,rp 6 C°°(SP). The functions in C°°(SP) have the following strong boundedness 
property 

(4) (rpd x y(r P d y y(r P d z ) k u £ C°°(EP) C X°°(P) 

for all u e C°°(SP). The consideration of C°°(SP) and of the derivatives of the 
form r P d x , rpd y , and rpd z is a substitute for the results on Lie manifolds used in 
PP. However, the results of pQ also apply to non-compact manifolds and to a larger 
class of singular domains. 

The methods of this paper are used for a general regularity and well-posedness 
result for anisotropic elasticity in general polyhedral domains (including cracks) in 
[4^j . We do not include in this paper any concrete applications, but let us refer the 
reader to El El CHI EH 122] , where concrete applications of results similar to ours 
were provided. 
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2. Smooth functions and differential operators on P 

In this section, we shall introduce the space C°°(SP) C C°°(¥)nL 00 (F) and relate 
it to the differentials i~rd x , rpd y , and rpd z mentioned in the Introduction. Similar 
vector fields have appeared also in ^7j. When only edges are involved [i.e., no 
vertices), the use of these vector fields goes back to ^2E2- See also |44l I50| . 

2.1. Polygons and polyhedral domains. Let us fix some terminology to be used 
in what follows. 

A polygon Po in a two dimensional Euclidean space is an open, connected subset 
whose boundary consists of finitely many straight segments (possibly of infinite 
length) called sides and having at most the end points in common. For simplicity, 
we assume that 8Pq = <9Po, which means that no point of the boundary dPo is in 
the interior of Po (thus cracks are excluded). The points common to more than 
one straight segment of the boundary are called the vertices of Po- We require that 
each vertex belongs to exactly two sides. 

We do not require the boundary of Po to be connected. For simplicity, in this 
paper we also assume that the sides are maximally extended, so that they are 
not contained in larger segments contained in the boundary. This assumption is 
however not essential. 

Similarly, a polyhedral domain P C I 3 is a connected, open subset whose bound- 
ary satisfies dP = dP = Uj=i Dj and: 

(i) each Dj is a polygon contained in an affine 2-dimcnsional subspace of M. 3 ; 

(ii) the sets Dj are disjoint; 

(iii) a side of Dj is a side of exactly one other Dk- 

The vertices of P are the vertices of the polygonal domains Dj . The edges of P 
are the sides of the polygonal domains Dj . Hence an edge belongs to exactly two 
faces of P. For each vertex P of P, we choose a small open ball Vp centered in P. 
We assume that the neighborhoods Vp are chosen to be disjoint. 

We stress that, in our convention, both the polygons and the polyhedra are open 
subsets. We do not require these sets to be bounded in general, although this 
assumption is needed for some of our results. 

2.2. Useful functions and other notation. Assume, for the definition of r e , 6 e , 
and <f>p te in this subsection, that P is convex. If P is not convex, then we slightly 
change the definitions of these functions such that the new functions retain their 
behaviour around e, but will become smooth everywhere in space except one p. 
The modified functions <fip. e and 9 e will then be defined and smooth on P. We 
postpone the technical construction of the modified functions 0p. e and 9 e for the 
Appendix, in order not to interrupt the flow of the presentation. (Let us stress, 
however, that none of our results requires the assumption that P be convex.) 

Let us first recall from the Introduction that we have denoted by pp(p) the 
distance from p to the vertex P of P. Also, recall that we have denoted by r e (p) 
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the distance from p to the line determined by the edge e of P and by 
(5) r P := JJr e x JJpp, where f e := p~^ p^r e for e=[AB]. 

e P 

In the above formula, the products are taken over all vertices P and all edges e 
of P. The notation e = [AB] means that e is the edge joining the vertices A and 
B. If e = [A, oo), that is, if e is a half-line, then r e := p~^r e . Finally, if e is infinite 
in both directions (i. e., for a dihedral angle), we let f e := r e . 

Choose for each edge e a plane V e containing one of the faces Dj of P such that 
e C Dj. If x is not on the line defined by e, we define 9 e to be the angle in a 
cylindrical coordinates system (r e ,9 e ,z) determined by the edge e and the plane 
V e - More precisely, let q € e be the foot of the perpendicular from p to e. Then 
9 e (p) is the angle between pq and V e - Similarly, for each vertex P and edge e 
adjacent to P, we define 4>p,e(p) to be the angle between the segment pP and the 
edge e (except for p = P, in which case 4>p,e{p) is not defined). 

If P is convex, then the functions 9 e and <j)p :e are defined and smooth on P (recall 
that P is an open subset). They will be part of the spherical coordinate system 
(pp,6 e ,(f)p te ) centered at P. For P non-convex, this property will be enjoyed by 
the modified functions 9 e and <fip, e introduced in the Appendix. All the following 
definitions and constructions below are the same in the case of a non-convex domain, 
but using the modified 9 and 4> variables. 

We shall denote by 9 — (9 ei , . . . , 9 £r ) the vector variable that puts together all 
the 9 e functions, for e ranging through the set of all edges {ei, . . . , e r }. Similarly, let 
{</>i, . . . , 4> p } list all the functions 4>p tC , for all vertices P and all edges e containing 
P we shall denote by <f> — (<pi, . . . , <f> p ) the vector variable that puts together all the 
4>p : e functions. We then introduce the space W k '°° (EP) as the space of functions 
u : P — > C of the form 

u(x,y,z) = f(x,y,z,0,<f>) = f(x,y,z,9 ei ,. . . , 9 er , <f> x , . . . , <f> p ), 

,/G W^°°(Px (0,2^) r x (0,ir) p ). 

Thus / above has k bounded weak derivatives. We let C°°(EP) := f] k VF fe '°°(EP). 
The point of this definition is that, for example, 9 e is a smooth function on P that 
is not in W k ^(P) for k > 1. On the other hand 9 e e C°°(EP), by definition. 

One can show as in ^ II Tj that there exists a canonical Riemannian manifold 
E(P) such that C°°(E(P)) = C°°(EP), so our notation is justified. The construction 
of a space with this property is not very intuitive. However, at this point, we do 
not assign any significance to EP, which should be regarded in this paper just as a 
symbol. (Let us mention however, that, had we used curved boundaries, then the 
desingularizations EDj of the faces would have been necessary. See |43].l 

2.3. Vector fields and C°°(EP). We now establish several technical properties of 
the functions in C°°(EP), especially in relation to the vector fields (differentials) 
rpd x , rpd y: and rpd z . 

Let us notice first that it follows right away from the definition that C°°(EP) is 
closed under addition and multiplication. 

Lemma 2.1. Let P be a vertex ofW, then p P e C°°(EP). Similarly, let e = [AB] 
be the edge o/P joining the vertices A and B, then r e := pT r e £ C°°(EP). In 
particular, r r := ]J e f e x ]J p p P 6 C°°(EP). 
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This is proved using polar coordinates. Assume P belongs to the edge e, then 
pp = (sin 4>p, e cos 6 e )~ l x, where this is defined (x stands for the first component 
variable). Similar formulas for pp in terms of y and z then combine, using a 
partition of unity on P 3 \ {P} with functions in C°°(£P), to define pp globally as 
an element in C°°(SP). 

Similarly, f e — pa sin(^ e , so f e /pA is "smooth" near A. The same argument, 
together with a partition of unity, shows that r e € C°°(£P). Our result then follows 
from the fact that C°°(SP) is closed under products, by definition. 

Lemma 2.2. Let flip) be the distance from p to the union of the edges o/P. Then 
there exists C > such that C^ 1 -d(p) < r P < Cd{p) for all peP. 

This lemma is proved using the homogeneity properties of the functions $ and 
rp close to the vertices and edges of P. Using a compactness argument, it is enough 
to prove that the ratio 7>/i? is bounded and bounded away from zero in the neigh- 
borhood of each point. This allows us to assume that P is cither a dihedral angle or 
an infinite cone. If P is the dihedral angle < 9 < a, with a fixed, then r P /d = 1. 
If P is a cone with center the origin, let a t be the dilation with center the origin 
and ratio t. Then rp(ctt(p)) — trp(p) and -d(at(p)) = td(p). This shows that the 
ratio JTp(p)/i9(p) depends only aa.p/\p\. Furthermore, rp is a continuous function 
on the compact set P l~l S*" -1 , and the lemma follows from this. 

Lemma 2.3. We have that the functions r e d x 9 e , r e d y 6 e , r e d z 6 e , ppd x (f>p^ e , 
ppd v 4>p, e , ppd z 4>p, e , d x r e , d y r e , d z r e , d x p P , d y p P , and d z p P are all in C°° (£P) . 

To prove this, let us notice first that we can use any linear system of coordinates 
(x, y, z). In particular, for each of the above calculations, we can assume that 
our cylindrical or spherical coordinate system is aligned to the coordinate system 
(x,y,z). Then the result is simply an exercise in the calculation of the partial 
derivatives of the cylindrical coordinates 9 and r and of the spherical coordinates 
4> and p. 

Corollary 2.4. We have d x r P , d y r P , d z r v 6 C°°(EP). 

Proof. Let us concentrate on d x . We use the product rule to compute the derivative 
of rp. A summand containing d x pp is in C°°(SP) by Lemma [2.31 Let e = [AB\. 
The other products are obtained by replacing f e := p^ pg r e with 

dxiTe) = p'a 1 Pp^dxire) - p~ A 1 d x (pA)f e - Pp?d x {p B )f e - 

The factors of p^ 1 and p^ 1 then cancel out in the product defining rp and all the 
remaining factors are in C°°(SP) by Lemma \'2. '6\ □ 

Proposition 2.5. If u G C°°(SP), then the functions rpd x u, r P d y u, and rrd z u 
are in C°° (EP). 

Proof. This follows from Lemma IO and r P £ r e C°°(EP) n j o P C 00 (EP). □ 

Let us denote by Diff™(P) the differential operators of order m on P linearly 
generated by differential operators of the form 

u(r r d) a := u{r P d x ) ai (r P d y ) a2 (r P d z ) a3 , \a\ := a x + a 2 + a 3 < m, u G C°°(EP). 

We agree that Diff™(P) := C oc (EP) and we shall denote Diffg°(P) := \J m Diff™(P). 
In case of edges (no vertices), similar algebras were considered also by Mazzeo, ^2 
l%2*] . Algebras more closely related to ours appear in To get more insight into 
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the structure of Diff ^° (P) , we shall need two simple calculations that we formalize 
in the following lemma, whose proof is based on the fact that djiy> G C°°(£P). 

Lemma 2.6. Let A G R and let dj and d k stand for either of d x , d y , or d z . Then 
rp X (r P dj)r$ - r v dj = Xdj{r P ) 6 C°°(EP), and 

[r P dj,rpdk} ■= {r P dj)(rpd k ) - {r v d k )(r v dj) = d 3 (r v )r v d k - d k (r v )r v d 3 G DiffJ(P). 

Then we have the following simple but basic result. 

Proposition 2.7. We have Diff £(P) Diff ™(P) C Diffo +Tn (P) and hence Diff^(P) 
is an algebra. 

Proof. We shall prove by induction on k + m that Diff § (P)Diff (P) C Diff g +rn (P) . 
Indeed, if k + to — 0, then k — m = and the statement is clearly true because 
C°°(EP) is closed under products. Let us assume then that k + m > 0. We need to 
show that u(r P d) a v(r P d)P G Diffg +m (P) ifw,u € C°°(SP) and \a\ := ai+a 2 + a 3 = 
k, \fj\ := Pi +/3 2 +/?3 = to, where a = (ai,a 2 ,a 3 ) and /3 = (fa, fo, fa). 
If rn = 0, then the relation 

u(r P d) a v = u(r P df [rrdj(v)] (r r d) a " 

for suitable a', a" with \a'\ + |a"| = k — 1, together with the induction hypothesis 
and with Proposition [213 shows that u(r P <9) Q u <E Diffg(P). 

Let now to be arbitrary. We shall proceed by a second induction on m. The same 
argument as in the paragraph above allows us to assume that v = 1. We can also 
assume that the monomial (rpd) a (r P d)^ is already ordered in the standard way. 
Then, using Lemma l2.6l we commute r P dj, the last derivative in (rpd) a , with rpd k , 
the first derivative in (rpd)@ . Induction on k+m for the terms containing dj(rp)rpd k 
and d k {rp)rpdj and induction on to or k + to for the term containing (rpd k )(rpdj) 
then complete the proof of the fact that Diff § (P) Diff ^(P) C Diff„ +m (P). □ 

The above proposition gives the following useful corollary. 

Corollary 2.8. If P is a differential operator of order to with smooth coefficients, 
then r™P G Diff™(P). 

Proof. It is enough to show that rfd a G Diff™(P) if d a = d^d^d? 3 with |a| = to. 
We shall again proceed by induction on to. The case m = 1 is obvious. Let dj be 
the first derivative in 9", so that d a — djd a . Then Lemma f2. 61 and Corollary 12.41 
give 

r^d Q - (rpdjXrp-W) = -(to - l)dj(r P ){r™- l d a ') G Diff J 1 " 1 (P). 

Then Proposition shows that Diff J(P)Diff™" 1 (P) C Diff™(P). This and the 
induction hypothesis allows us to complete the proof. □ 

The proof of the above corollary also shows that 
(6) r^d^d^d^ - (rpd x ) a Hrrd y r 2 (r P d z ) a3 G Diff^P), \a\ = to. 
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3. Function spaces on P 

We now recall and study the Babuska-Kondratiev spaces /C™(P) := /C"'(P;??) 
and /C™(<9P; $) on a 3-simensional polyhedral domain P and its boundary <9P. These 
spaces are weighted Sobolev spaces with weight given by the distance to the set 
of edges of P, as in Equation (J2J). Note that we can replace $ with r>, by Lemma 
12.21 (we shall use this below). 

3.1. The Babuska Kondratiev spaces. We let 

(7) Wg'P' a (P) = {u : P -> C, rj, o| ~ o 0°u G L P (P), for all \a\ < k}, 

for k G Z+, a G M, p G [l,oo]. If p = 2, we denote £*(P) := W|£' a (P), which 
coincides with the definition in the Introduction (Equation |2J). 
We similarly define 

Wgg' a {dP) = {u:dP^C, 4- a P{u\ D] ) e L p (L>j), for all k < m 

and all differential operators P of order k on Dj, k < m}, m G Z + . 

We let £*(aP;#) := W^?' a (<9P). Thus /C*(<9P;tf) ~ ©fcjp^i?) is thus a direct 
sum of weighted Sobolev spaces. Note that we require no compatibility conditions 
for the resulting functions on the faces Dj. 

Equation © and Lemma 12 . 21 then give immediately the following lemma. 

Lemma 3.1. We have /C™(P) = {it, r fl Pu G L 2 (P), for all P G Diffg(P)}. A 
similar result holds for /C™(<9P; i?) and /or W|'£' (P). 

Next, Proposition 12. 51 and Corollary 12.41 together with a straightforward calcu- 
lation, show the following. 

Lemma 3.2. The multiplication map Wg^' b xK™{P) 3 (u, f) i-> uf G £™ +fe (P) * s 
continuous. We also have C°°(EP) C VF^°°'°(P) and r P G VF^°°' b (P), and Zience 
the map /C™(P) 9 u i— > r$,tt G A^|_ b (P) is a continuous isomorphism of Banach 
spaces. 

From this lemma we obtain right away the following result. 

Proposition 3.3. Let k > m. Each Pq G Diff™(P) defines a continuous map 
Pa • Kj(P) — ► /C^ _m (P). XTie family r^^Pgrp is a family of bounded operators 
/C„(P) — > A^ _m (P) depending continuously on A. 

Similarly, if P is a differential operator with smooth coefficients on P, then 
Tp A Prp defines a continuous family of bounded operators /C„(P) — > /C^Z™(P). 

Proof. The first part follows from Lemma |3.1I The second part follows from the 
first part of this proposition and Lemma 13.21 □ 

o 

We define the spaces /C~ fc (P), k G Z+, by duality. More precisely, let /C„ (P) be 

o 

the closure of C C °°(P) in £*(P). Then we define £l*(P) to be the dual of (P), the 
duality pairing being an extension of the bilinear form (u, v) i— * J p uv dvol. With 
this definition, we can drop the requirement that k > m in Proposition 13. 31 

Let us also note that the resulting weighted Sobolev spaces on the polygons 
Dj are different from the weighted Sobolev spaces obtained by using the distance 
to the vertices of these polygons. A regularity theorem on Dj would involve the 
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latter weight (as in Kondratiev's theorem mentioned in the Introduction). A 
consequence of this is that the spaces K.\ (dP; &) behave more like the Sobolev 
spaces defined on a smooth manifold without boundary than like the Sobolev spaces 
defined on a bounded domain with (smooth) boundary. In particular, we define 
KT£(<9P;i?) as the dual of K.%(dT>;<&). The spaces /C*(<9P), s g Z, can be defined 
by interpolation, although in this paper we shall use a different definition using 
partitions of unity (see the following subsection; the two definitions are equivalent, 
although we shall not need a proof of this fact in this article). 

3.2. Definition of Sobolev spaces using partitions of unity. As in pQ, it is 

important to define the spaces A^(P) using partitions of unity. Similar construc- 
tions were used in j^nHSES ED This construction is possible because the spaces 
/Qw 2 (P) are the Sobolev spaces associated to the metric r^ 2 gE, where gE is the 
usual Euclidean metric. 

We shall need the following lemma. Recall that d(p) denotes the distance from p 
to the edges of P. In view of Lcmma l2~2l in all estimates involving i9, we can replace 
$ with rp, although not the other way around, because i9 is not smooth. 

Let <9 sing P be the union of the edges of P and P' := P \ <9 sing P. 

Lemma 3.4. There is eo G (0, 1), an integer k, and a sequence C m > of constants 
such that, for any e G (0, eo], there is a sequence of points {xj} C P' := P \ d S mgP 
and a partition of unity (f>j G C£°(P') with the following properties: 

(i) either B(xj,ed(xj)/4) is contained in P or Xj G dP, $(xj) > 0, and the ball 
B{xj,e${xj)) intersects only the face Di to which Xj belongs; 

(ii) supp(0j) C B(xj,e-d(xj)/2) if xj G dP and supp(0j) C B(xj,efl(xj)/8) oth- 
erwise; 

(iii) 4>j{xj) = 1 and ||(rp9)°^-||i«»(p) < C\ a \e^ ; and 

(iv) a point igP can belong to at most k of the sets B{xj, e$(xj)). 

Let us notice that B(xj, ed(xj)) does not intersect any edge of P because e < 1. 
Moreover, the conditions that ||r>V<?!>j||z,oo < C\eT x and 4>j(xj) — 1 guarantee that 
the support of 4>j is comparable in size with ed{xj). This is reminiscent of the 
conditions appearing in the definition of the Generalized Finite Element spaces 

izi Ei nni - 

A proof of this lemma will be given in the Appendix. It is essentially a result 
that, in the case of non-compact manifolds, goes back to Aubin. It was subsequently 
used by Gromov and in 1511 1521 15*%| . We shall fix e = eo m what follows and a 
sequence Xj and a partition of unity (j>j as in the lemma. 

Lemma 3.5. Let Uk — Ylj=i 4>j u > f or u £ ^-^(P)- Then — » u in /C™(P). 

Proof. Let <&k '■= Si=i We have that the sequence (rp9)"$fe is bounded in 
the 'sup'-norm and converges to pointwise everywhere if a ^ 0. Similarly, <frfc is 
bounded and converges to 1 pointwise everywhere. The result then follows from 
this using also the Lebesgue dominated convergence theorem. □ 

Denote by ctj(x) = Xj + $(xj)(x — Xj) be the dilation of center Xj and ratio 
$(xj). Let J be the set of indices j such that Xj G dP. Below, by H m we shall 
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mean either H m (R 3 ) or H m (R%). Also, denote by 

(8) v m , a {u? ■■=J2^) 3 ~ 2a W^ u ') oa ^ 

j 

■= £^) 3 - 2a ii(<M ° M W 3 ) + E^) 3 ~ 2a n(<M ° M W^r 

We agree that ||(0/it) ° ctj\\ H ,n = oo if (4>ju) o ctj £ H m (M. 3 ) (or if (<f)ju) o ctj E - 
H m (M. 3 ,), respectively). Note that the functions {(f>ju) o atj will all have support 
contained in a fixed ball, namely, the ball B(0, e /2) of radius e /2 and center the 
origin. Moreover, all derivatives d a (4>j o ctj) are bounded for each fixed a and 
arbitrary j by Lemma 13.41 

Proposition 3.6. We have u E /C™(P), to £ Z, «/, and on/y z/ ; v m ,a(u) < oo. 
Moreover, v m>a (u) defines an equivalent norm on /C™(P). 

The proof of this Proposition is standard (see |KSI Lemma 2.4], £Q, or |52]); for 
m < one also has to check that both definitions are compatible with duality. We 
include a brief sketch below. 

Proof. Let us also introduce 

V m ,a{ U f ■= X! W(^ U )W^(Fy 
3 

Then the fact that d(x)/-d(xj) and d(xj)/d(x) are bounded by (1 — e) _1 on the 
ball B(xj,ed(xj)), for e E (0,1) and a change of variables shows that v m>a and 
Vrn.a define equivalent norms. It is then enough to prove that v m ,a(u) defines an 
equivalent norm on /C™(P). For m = 0, this follows from the inequalities 

\\ r r a u\\ 2 L 2 {P) < Kvl a < n\\rp a u\\ 2 L2m . 

For arbitrary to, we use induction on to and the fact that Y]j \ (rpd) a 4>j(p)\ is 
bounded uniformly in p E P for all a. □ 

We proceed in the same way to study the spaces ZC*((9P;t9), s E R. Let us 
identify the plane containing each face Dk of P with a copy of R 2 . Then let 

(9) p s , a {u) 2 :=^tf(^) 2 - 2 l(<M°a;llW 2 )' seM +' 

Note that only the indices j for which Xj E 9P are used above. Also, note that 
the power of $Xj was changed from 3 — 2a to 2 — 2a. 

Then we have an analogous description of the spaces K.^(d¥;-d), s E Z. 

Proposition 3.7. We /iawe u E /C*(<9P; #) if, and only if, /i s a (u) < oo. Moreover, 
/j, Si a(u) defines an equivalent norm on /C*(9P;??), s E Z. 

We can then define /C*(<9P; ■#), s E R, as the space of functions u for which 
A t s,a(u) < oo with the induced norm. From this we obtain, by reducing to the 
Euclidean case, the following Trace Theorem. Let <9 s ingP be the union of the edges 
of P and P' := P \ <9 sing P, as above. 

Theorem 3.8 (Trace theorem). The space C~(P') is dense m £™(P), m E Z+. 
The restriction to the boundary extends to a continuous, surjective map /C™(P) — > 
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£™_i/2 f or 771 > !■ -for m = 1, ifte kernel of this map is the closure of 

C~(P) w/Cj(P). 

Proof. Clearly C C °°(P') C /C™(P), for any m £ Z+ and any a 6 R. To prove that 
it is a dense subspace, let u G /C™(P). By Lemma [3.51 we may assume that the 
support of u does not intersect <9 s i ng P (replace u with for some k large). Then 
we use the fact that C°°(fl) is dense in H m (Q) for Q, a smooth, bounded domain 
and the fact that the H m -norm is equivalent to the norm on /C™ (P) when restricted 
to functions with support in a fixed compact K such that K does not intersect any 
edge of P (i. e., K n <9 sing P = 0). 
We have 

Mm-l/2,a-l/2(w|ap) 2 := ^tiiXj'f^Wi^jU) °Q!j|ap||^m-l/2 (K 2) 

< C^^(^) 3 - 2a ||(^ u )oa,||^ m(R3) < Cu m , a (u), 
je.J 

and hence the restriction map /C™(P) -» K™^* {&¥]<&) is defined and continuous 
for m > 1, by Propositions 13 . 6l and 13 . 7l To prove that this map is continuous, let us 
fix a continuous extension operator E : ii" m_1 / 2 (R 2 ) — > iP"(R 3 ). By rotation and 
translation, we extend this definition to an extension operator E : 7J m_1 / 2 (y) — » 
ff™(l 3 ), for any two dimensional subspace Fcl 3 . 

Let then v : <9P — > C be a function in KJ^_S^{d¥]'&). Let us fix a function 
ijj G C(?°(R 3 ) with support in the ball B(0, eo) of radius eo and center at the origin 
such that ip = 1 on P(0,eo/2). Let vj{p) = (f>j(oij(p))u(otj(p)), which is defined on 
a subspace of R 3 of dimension 2. We define 

3 

Then u G /C™(P) and u\ 9r = v. 

Finally, let u G /C*(P) such that = 0. Let u& be as in Lemma 1531 Then 
Uk\dP = 0. Using again the equivalence of the H 1 and /C„(P)-norms on functions 
with support in a fixed compact set K such that if n9 s i ng P = 0, we see that each Uk 
can be approximately in /C*(P) as well as we want by a function Vf. G C£°(P'). Then 
we can take our approximation of u to be v = Ylk=i Wfc > ^ or ^ large enough. □ 

4. Proof of the regularity theorem 

We include in this section the proof of Theorem ll.il Its proof is reduced to the 
Euclidean case using a partition of unity 4>j satisfying the conditions of Lemma 13.41 
for e = eo, as in the previous section. 

Proof, (of Thcorcm ll.lO The trace theorem, Theorem l3 . 81 allows us to assume that 
u\dP = 0. We then notice that, locally, Theorem 11.11 is a well known statement. 
Namely, let us consider a function v with support in the ball of radius eo- We 
assume that either v G iJ 1 (R 3 ) or v G Hq(M.+ ) (that is, v = on R 2 , the boundary 
of R 3 ^ = {z > 0}). Then there exists a constant C > such that, for all m > 0, 

(io) hf Hm+ , m < a(||A«||^-i (B 3) + IMli^aj). 
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or, respectively, 

(ii) INI <C r (\\Av\\ 2 Hm ^ (w3+) + \\v\\ 2 L2{R3+) ). 

The constant C r in the two equations above depends only on eo. (In fact, Equation 
(|10|l implies Equation |H) , by taking v to be odd with respect to the reflection in 
the boundary of the half space R+.) 

We shall proceed by induction on m > 0. For m — 0, the result is tautologically 
true, because of the term ||«||je° x (p) on the right hand side of the regularity esti- 
mate of Theorem ll.il Let now {<fij} be the partition of unity and aj be dilations 
appearing in Equation ©. In particular, the partition of unity <j)j satisfies the con- 
ditions of Lemma f3. 41 which implies that supp(^j) C B(xj , eo'd(xj)/2) if Xj G <9P 
and supp(^ij) C B(xj,eo'&(xj)/8) otherwise. We also have that all derivatives of 
order < k of the functions <fij o ay are bounded. This implies in turn that the 
commutator 

Pj := [A, o aj ] := A(<j>j o aj ) - {4>j ° aj)A 

is a differential operator all of whose coefficients have bounded derivatives. 

Let |H|ff™ denote either |H|,ffm(K 3 ) or ||«||ijm(g3 ), depending on where the func- 
tion v is defined. Let r\j — ip o aj 1 , where ip € C^°(IR 3 ) has support in -6(0, eo) and 
is equal to 1 on J5(0, eo/2), as before. 

Then Equations l|10[l and (|11|1 and the above remarks give 

fm+2,a(w) 2 := ^tiixrf-^W^jU) O Ujfjj^ 
3 

< c r J2^) 3 - 2a (\\M(M) °^]f Hm + IK^o^iil.) 

3 

< Cj2^jf- 2a (\\(<l>joaj)A(uoaj)\\ 2 Hm + \\Pj(u o aj)\\ 2 Hm + ||(<^) o aj\\\.) 

3 

< C^ 7 9(^) 3 - 2a ( 7 9(^) 4 ||(^A^)oo,)||| fm + ||( %u )o aj -||^ m+1 + 11(0^)0^111,) 

3 

< C(v m ,a-2(Au) 2 + ^ V m +l,a(V3u) 2 + V , a (u) 2 ). 

3 

Since no more than k of the functions rjjU axe non-zero at any given point of P 
and all the derivatives (rpd) a T]j are bounded for all fixed |a|, we obtain that 
J2 3 v m+ia(i),i) 2 < Cis m+1 , a (u) 2 . This then gives 

(u) 2 < C (ls m , a - 2 (Au) 2 + V m + lA u ) 2 )- 

By induction on m we then obtain 

V m +2,a{ U ) 2 < C{u m . a -2{Au) 2 + h> ,a{u) 2 ) . 

The result then follows from Proposition ^. 61 which states that the norms || • ||x;* rp) 
and vt, a are equivalent. □ 

See ^2] f° r applications of these results, especially of the above theorem. 
By contrast, it is known that in the framework of the usual Sobolev spaces 
# m (P), the smoothness of the solution of Q is limited fZZ\ E3 EE1 E01 E3 
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Appendix A. Additional constructions 

In this appendix we explain how to modify the constructions of the functions 
9 e and <pp ie introduced in Section [21 when P is not convex and how to construct a 
partition of unity satisfying the conditions of Lemma 13.41 

A.l. The modified functions 9, cf>, and r e . We continue to denote by pp(p) the 
distance from p to the vertex P. By a dilation, we can assume that each edge of P 
has length at least 4. 

Let us first modify the functions cf>A,e- We can find S > small enough so that 
for any vertex P, the sets 4>P.e > tt — 2<5 do not intersect (e ranges through the set 
of edges containing P). Let e = [AB] and tpi : [0, tt] — ► [0, 1 — 6] be a smooth, non- 
decreasing function such that ipi( x ) = x f° r < x < tt — 26 and ipi( x ) = tt — <5 for 
x > tt — 6. Also, let ip2 '■ [0, oo) — > [0, 1] be a smooth, non-increasing function such 
that ip2(x) = 1 for < x < 1 and ip2( x ) — for 2 < x. Then we replace 4>A,e with 
, 0i(0A,e)'02(p J 4)- This modifies the function (j)A,e to make it smooth everywhere 
except on e. 

We now modify the functions 9 e . They will be modified in two ways. Let us 
fix an edge e = [AB]. To understand these modifications, it is useful to think of 
the spherical domain loa associated to the vertex A. The old function 9 e served 
the purpose of both desingularizing loa close to the vertex associated to e and of 
providing global coordinates on loa away from the vertices (toghether with the 
functions 4>A,e)- These two purposes of the old 9 e will be accomplished by two 
modified functions 6. Let ^3 : [0, tt] — > [0, 1] be a smooth, non-increasing function 
such that ip3(x) = 1 for x € [0,a] and ip3(x) = for x > 2a. We then similarly 
modify 9 e by replacing it with ip3(4>A,e)' t l ) 3(4'B,e)' t p2{r e )0 e . This will make 9 e defined 
and smooth everywhere in space except on e (if 7 is large enough). The resulting 
function 9 e serves the purpose of desingularizing u>a near the vertex corresponding 
to A. Let next ip4 : [0, 2tt] — ► [0, 2w] be a smooth function such that ip4 (t) = t for 
t e [2e, 1 - 2e] and xp 4 (t) = for t G [0, e] U [1 - e, 1] . The second kind of functions 9 A 
will be obtained by considering V , 4(#e)V'4(20 J 4,e)V ; 4(20,B,e) for e > small enough 
and all choices of faces passing through e. (To define the old functions 9 e , we first 
chose a plane through e and containing one of the faces of P. This plane was the 
plane where 9 e = 0. For the new functions 9 e , we consider all the planes through e 
and containing one of the faces of P.) These new functions will be smooth on u>a 
near its vertices, but provide global coordinates away from the vertices. 

Finally, let ^3 and e = [AB] be as in the above paragraph. We then replace r e 

with 1p3,{(j>A,e)llj?,{(t>B,e)re + (1 - ^3 {4>A,e ))PA + (1 ~ 4>3 {<t>B,e))PB ■ 

Let us notice that one can define SP directly, which would provide the definition 
of C°°(SP) as the space of smooth functions on SP pQ. The advantage of the 
approach in |J is that it makes no distinction between the cases when P is convex 
or non-convex. The approach in this paper has the advantage that it is much 
simpler and more intuitive in the convex case. 

Further intuition in the construction of the spaces C°°(5]P) can be obtained 
from the paper |22) . page 254, by Costabel and Dauge where various regions and 
subregions of a polyhedral domain were analyzed. See also 1191 l3"fi| . 

A. 2. The partition of unity. Our partition of unity will depend on parameters 
(a, b, c) that will be specified below. 
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First of all, let us denote by B(P, n) the open ball of center P and radius 2~ n a. 
By choosing a small enough, we can assume that the balls B(P, 1) do not intersect. 
Then let E e>n be the set of points pGP that do not belong to any B(P, 2) and are 
at distance < 2~ n ab to the edge e. By choosing b small enough, we can assume 
that the sets B e ,i do not intersect. Let fli be obtained from P by removing the sets 
B{P,2) and E e #. 

For each edge e, let N e be the plane normal to e. Project E e> i \ E e ^ onto N e . 
The projection will be the intersection of an annulus with an angle. Denote this 
projection by C e . We shall cover C e with disks of radius c/2 and with disks of 
radius c/8. The disks of radius c/2 have the center on the straight sides of C e (the 
ones obtained from the angle) and the disks of radius c/8 that have centers in the 
interior of C e at distance at least c/4 to the angle defining C e . This yields the disks 
Di,..., Dn with centers q\, . . . , qtj. 

Let z be the variable along the line containing e. Then we cover E e ,k+i \ E e ^+2, 
k G Z + , with balls of radius 2~ k c and centers of the form (2~ k qj, 2~' £ ~ 3 c), if 2~ k qj 
is on one of the faces of P and is inside E e \. Otherwise, we consider the ball of 
radius 2~ fc ~ 2 c with centers of the form (2 qj, 2~ k ~ 3 c) as long as the center is still 
inside E e% \. 

Let us cover 

fii :=P\ (|jB(P,2)u|J J B e)2 ) 

P e 

with finitely many balls of radius c/2 or radius c/8 with centers in fij such that 
the balls of radius c/2 have the centers on the faces of P and the balls of radius c/8 
are at distance at least c/4 to the faces of P. 

Let Dp i, . . . , Dp : M : ... be the balls already constructed with centers in B(P, 1) \ 
B(P, 2). Then consider also the balls 2 Dpi, . . . , 2~ fe £>p ! jv, ■ ■ ■ obtained by dila- 
tions of ratio 2~ k and center P. We repeat this construction for all vertices P and 
all k G Z_|_. We consider all the balls D±, D2, ■ ■ ■ , constructed so far (relabeled into 
a sequence) from the coverings of E e< i, Oi, and from the dilations of ratio 2~ k for 
all the vertices P, as already explained. If we choose c small enough (after the 
choices of a and b have been made as explained above) , then the sequence of these 
balls is locally finite, the center of each ball is either on the faces of P or the closure 
of the ball is inside P. Moreover, for any such ball D with center p and radius r, we 
have that r/-&(p) is bounded from above and bounded from below from zero, say 
r/d(p) G [e , £q 1 ], for some eo G (0, 1). There is an integer k such that no k + 1 of 
the balls constructed have a common point. 

To any ball D of center q and radius r we associate the bump function 4>d{p) '■= 
ip(\p — q\/r), where ip : [0, 00) — > [0, 1] is smooth, is equal to 1 in a neighborhood 
of 0, is equal to in a neighborhood of [1, 00), and is > on [0, 1). Then we let 
T] = X^V'-Dj an( i 4>j — "0-Dj/ 7 ?- By further decreasing c, if necessary, we see that our 
partition of unity (together with the points Xj obtained as the centers of our balls) 
satisfies the conditions of Lemma 13.41 for the eo chosen above. 
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